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Abstract 



Using a regularisation of the chiral SU(2) Yang-Mills model by an 
C^il infinite number of Pauli-Villars fields the doubler free gauge invariant 

lattice model is constructed. Its continuum limit provides the recently 
proposed mechanism for cancelation of Witten anomaly. 



^ 1 Formulation a model containing chiral fermions on a lattice is a long 

standing problem due to the phenomenon of fermion doubling It was 
proposed by Wilson to cure the doubling by addition of an extra term to the 
naive lattice action which produces large masses of the order of cut off to 
doubling states. This eliminate doublers in continuum limit but it destroys 
the chiral invariance as well. It also is the source of chiral anomaly in the 
continuum limit. In some cases, however, such as QCD in the limit of zero 
mass u and d quarks, SM or other nonanomalous models this destruction 
is too "strong". In this case one still needs additional "fine tuning" of the 
model in order to approach the chiral limit. The counter-terms needed for the 
"fine tuning" are not known a priori, and require a careful (non-perturbative) 
computation (for a recent review see e.g. 0). From the other hand, no local 
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lattice formulation preserving the non-anomalous chiral symmetry for finite 
lattice spacings is known in present. 

Moreover, there exists a theorem due to Nielsen and Ninomiya [|3] which 
states that under certain circumstances such as locality, hypercubic and 
gauge invariance and hermicity there are no chiral fermions on a lattice. 

From the other hand among the nonanomalous models there exist ones 
given by the real representations of gauge groups. For such models one can 
formulate a lattice description satisfying all conditions of Nielsen-Ninomiya 
theorem but with no doubling in continuum limit. This is possible because 
both chiralities in real representations are equivalent. For example for a 
model describing a chiral fermion in n-dimensional representation of SO(n) 
one can write down the gauge invariant Wilson term as follows 

-aV> T CAV> + h.c, 

where a is the lattice spacing, C is the Dirac charge conjugation matrix, A 
is the lattice Laplace operator and "T" stands for transposed spinor. 
In case of 577(2) model described by the Lagrangian 

C = ~F ia , + 4+V1>+ (1) 
P_V+ = 0, P±=i(l± 7B ) 

where ip + belongs to the fundamental representation of SU(2) one can con- 
struct a field 

^ = V?Cz<7 2 , (2) 

where C is the same Dirac charge conjugation matrix and o"2 is Pauli o^- 
matrix (here it plays the role of 577(2) charge conjugation matrix). One can 
see that if), defined in such a manner has the same chirality but transforms 
by conjugate action of the gauge group as initial field One can try to 
write down the gauge invariant Wilsonian term in analogy with SO(n) case 
but since the product C ■ 102 is symmetric this term will be identically zero 
due to anticommuting properties of the fermionic field. As we will see later 
it is connected with Witten anomaly. If the number of fermionic fields were 
even there would be no problem since one can take the crossing product of 
different fields solving the problem. But one can re-define the fields in this 
case to show that in fact one deals with half number of Dirac fermions. 
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In what follows we will adopt a different strategy. In continuum one can 
regularize this model by introducing an infinite number of Pauli-Villars (PV) 
fields H as follows 

oo 

= ^+VV>+ + J3? r (V + AfrM- ( 3 ) 

r=l 

where ip T are Dirac spinors with Grassmannian parities (— ) r . Regularizations 
of this type were used to prove the gauge invariance of the continuum limit of 
lattice models with Wilson fermions 0, Smit-Swift model || and for chiral 
gauge invariant lattice formulation of the SO (10) unified model with SLAC 
fermions ]7|Q. The same regularization was used recently for the construction 
of the representation for global Witten anomaly in continuous SU(2) chiral 
YM model 0. 

Let us return to the lattice model. We will introduce the following regu- 
larized lattice action 

£re, = ? + ( / B + ^j{5,A} N )«) + + (4) 



£{^(B + ^{B,A} 

r=l k v 



tj) r + Mrip r ip r 



where D is the naive lattice Dirac operator given by 

(x) = i J2 ^ (^WO + ae ») - i> (*)) ( 5 ) 



here a is the lattice spacing and is the unite vector in the x^ direction, 
U^(x) is the lattice gauge field corresponding to the link connecting x and 
x + ae^, A in the eq.(f|) is the lattice Laplace operator the same (up to an 
extra factor a -1 ) as used in the Wilson term and M and A are new cut off 
parameters. The action (M generates the following Feynman rules 



P 2 (p) + M 2 r 2 : 



r^O (6) 



alternative way may consist in considering second scale regularization using a 
coarser lattice, for a recent review see M, and references therein. 
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for propagator and 

for the vertices. In the above equations P a (p) is given by 

Pa{p) = ~ sinp^a ( 1 + — — V(l - cosp Q a) ] . (9) 



a 



If n — 1 in the vertex we call it a simple one and if n > 1 we call it multiple 
one. 

Let us show that in the framework of the perturbation theory one will 
have the contribution of doubling states vanishing in the continuum limit if 
the cutoffs M and A go to infinity as 4^. 

To analyze the continuum limit let us show as an example calculation 
of a diagram with n simple external legs. Generalization to multiple legs is 
straightforward. The Feynman integral corresponding to the diagram looks 
as follows 

i%;X(ki,..,K) = r^Jl -^trP+v^p + kMp + h)... (io) 

a 

oo 

^(-) r tr1/ Ml (p,p + h)s r (p + kx)....V^{p - k nj p)s r (p) 

r=l 

where T ll ~- ln is the group factor. 

Function P a {p) given by the eq.(§) near the edge of the Brillouin zone 
behaves as follows 

Pa(p')^^p' a + 0(1), (11) 

where 2£ = + (— ) nfl ), and quantities n M and p' are defined by the 

following 

71 . 71 71 

2a Ffl ^ a 2a 
At the center of the Brillouin zone P a {p) has the usual expansion 

P a (p) ~Pa + •■• 
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where dots stay for terms vanishing as a 2 and faster. 

The propagator and the vertex at the edge of the Brillouin zone behaves 
like this 

S ' (p,) *> -SjhWr ' P, " P '" 7 " (12) 

One can see that if there were no PV regularization one would have near 
the edge of the Brillouin zone almost the same integrand (the factors 

A 2 2 

cancel out) except the mass term which here vanishes as -j^-M — > 0. 

Since one have the PV regularization the situation changes dramatically. 
Indeed, the leading part of the integrand behaves like (hereafter we omit 
primes for pf), 

~ BH > (p) ; M , r2) j ^))^>)~ 

— — I (P'(p)) n p n (p) 

dP 2 {p)J MP 2 {p)smh{7iP{p)/M) y KPn {P) 

The regularization produces exponential cut off for zones where P(p) 3> M ~ 
As lattice spacing goes to zero and A ~ M ~ ^ to infinity one will have 
the cut off "starting" from values of p ~ and "finishing" by the values 
of p approaching the Brillouin zone edge as \p — - \ ~ y/a. Indeed from eqs 
( fnirfj ) one can see that for lattice spacings a small enough the asymptotic 
behavior of the integrand near the Brillouin zone edge is as follows 



2£ dp 2 ] M p 2 smh(2^7r 2 p/MA 2 a 2 ) 
while in the center of the zone integrand behaves like 



(15) 



d\ 71—1 
\ 71 



dp 2 J Mp 2 sinh {n 2 p/M) 



(16) 



From the above equations one can see that the integral over central zone 
tends to the correct continuum value (after extraction of UV divergencies) 
while integration over the strip close to the edge of the Brillouin zone decrease 



5 



with a going to zero as ~ a 2 . In the domain of momenta between these two 
extrema i.e. when p is departed from the center by more then value of M 
and from the edge by MA 2 a 2 ~ y/a the integrand decays exponentially and 
its contribution to the integral vanishes in the continuum limit. 

In fact to be rigorous in the analysis of the behavior of the Feynman 
integral near the edge of the Brillouin zone one cannot neglect the exter- 
nal momenta k since the loop momentum p here becomes small. But as 
can be shown including external momenta in the analysis cannot affect the 
conclusion. 

As a result for any diagram in perturbation theory one have the contribu- 
tion of the doubling states vanishing in the continuum limit. 

From the other hand it is known that SU (2) chiral YM model suffers 
from a global anomaly ||10|| . As it was shown by Witten one cannot globally 



define the sign of the fermionic determinant once topologically nontrivial 
gauge transformations which can change the sign of the determinant are 
considered. 

In our case, however, it seems everything to be all right. One has gauge 
invariant regularized action with doubler's contribution suppressed in the 
continuum limit. And no gauge noninvariant effect like change of the deter- 
minant sign can be produced here. The problem is that the above constructed 
lattice regular izat ion is equivalent to the model considered by Witten only 
perturbatively. 

To illustrate this consider eigenvalue problem for (modified) chiral Dirac 
operator V = (^D + ^2 {D, A} J . Charge conjugation relates both chiralities 
in such a way that one can introduce scalar product such that chiral Dirac 
operator is symmetric. Let {Aj} be the set of its eigenvalues. Then one has 
for the following expression for the determinant (see ref.||) 

lel *=Ilx-(») (17) 

one can see from this that the large eigenvalues (|Aj| 3> M) are cut off but 
not their signs. The same happens with the contribution of the doubling 
states. Arguments from the perturbation theory shall tell us that J^j — > 00 
where Aj are eigenvalues corresponding to unwanted (doubling) states. But 
when they are cut off they leave their signs alive in eq. fll7D . From the other 
hand from general reasons (e.g. Nielson-Ninomiya theorem) we know that 
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total number of the minus signs in eq. flTT] ) is an even oneQ. This would 
prevent the determinant from being negative. 

Indeed, for a sufficiently small lattice spacing one can write down the 
regularized determinant as follows (see Appendix) 



detP « ff^tanh ^ X (18) 



x 



JJ JJjWI t ail h 



, , , , , ^ V 2M* 

{doublers} J 

where are eigenvalues of the continuum (chiral) Dirac operator and M* 
is some effective cutoff mass of the order MA 2 a 2 ~ y/a. We wrote the 
regularized Dirac operator determinant as a product of determinants of 16 
Dirac operators. Each of them goes to ±1 as a goes to but their product 
is always +1. 

This lead to vanishing of the Witten anomaly as it is seen from the lattice. 
In fact on a lattice there is no reason for Witten anomaly since the topology 
of the lattice is different from one of the continuum. 

Now let us note that the cancelation of the phase factor here is identic 
to the mechanism proposed by Slavnov in ref . but in this case cancelation 
arises "naturally" from the lattice. The role of fields introduced in ref. to 
cancel Witten anomaly is played here by the lattice doubling states. 



Discussion 

In the present work we proposed a chiral gauge invariant lattice regulariza- 
tion of SU(2) YM model with correct (perturbative) continuum limit. The 
method includes using a second scale regularization with infinite number 
of Pauli-Villars fields |3| and insertion of a high order derivative term into 
fermionic action. 

There is no doubling or other lattice artefact on the perturbative level. 
But despite vanishing perturbatively, the doubler's contribution produces the 
Witten anomaly cancelation mechanism proposed in ||. 

Another lecture we learned from this is that a lattice model which is 
perturbatively equivalent to some continuum model can, however, be dif- 
ferent from it at nonperturbative level. Moreover, since the only way of 

2 Moreover, due to hypercubic invariance it is oc 2 4 = 16. 
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defining nongaussian path integral, beyond the perturbation theory is the 
lattice discretization one is tempted to define nonperturbative features of 
the continuum model as an extrapolation of the lattice one. 

Since for a given lattice spacing the momentum function P M (p) is bounded 
from above one can limit oneself to a finite number of PV fields since the 
fields with Mr 3> - will decouple. In this case decay of the integrand near 
the edge of the Brillouin zone will be polynomial rather than exponential. 
The same arguments and conclusions, however, will remain valid for this 
choice. But for practical calculation it seems to be more convenient to use 
the Grassmannian even path integral where the gauge invariant Wilson term 
can be constructed for calculation of inverse (absolute value of) fermionic 
determinant. 
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Appendix 



In this appendix we will derive the expression (18) for the determinant. Since 
the gauge fields are external ones for sufficiently small lattice spacings 
a they can be chosen close to the unity. Also let T A , (A = 1, ...,15) be 
generators of the symmetry which relates fermionic doubling states (see e.g. 
TO 



{T A } = {I, T M , FT" {y>y), t x T 2 T 3 T 4 }, (Al) 

where 

T" = W s(-)^, 
commutes with the naive Dirac operator 

DT A = T A D 

This symmetry is responsible for fermion doubling. Laplace operator -hA is 
not invariant under this symmetry^ but it transforms as follows 

T~ l AT a = A + {U a {x)ijj x+aSl + £/_ a (x)^-as) (A2) 



This is also the reason why Wilsonian term kills fermion doubling states. 
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the second term under our conditions is approximately j^2?p x . 

Now consider the modified Dirac operator D = (^D + -^{D, A}j and 

consider eigenvalues of T> lying in the subspace of small eigenvalues of A, 
A(A) <C 1. The spectrum of the modified Dirac operator T> here is close to the 
spectrum of the naive one (and for small lattice spacings also to the spectrum 
of the continuum Dirac operator). If now to act on such an eigenstate with 
eigenvalue fii by a transformation T the resulting state will have eigenvalue 
close to jr^fii- This way one finds fifteen such states. For greater values of 
D this relation becomes less and less exact but greater values are cut off and 



they do not contribute (except the sign) to the determinant (|T8|) . The extra 
factor ~ j^2 a 2 can be included in the cut off mass by its rescaling 

M -> M* = MA 2 a 2 . (A3) 

Finally let us note that the requirement that — > 1 as a — > is impor- 
tant here. This is natural since gauge fields are external. But if one wants 
to consider dynamical gauge fields here one has to introduce a regulariza- 
tion with a cut off of the order ~ \j\fa also for the gauge fields. This can 
be easily seen from the perturbation theory if one considers e.g. self-energy 
Feynman diagram for the fermionic field. Since PV mechanism cannot be 
applied here one has to have a cutoff for the gauge fields to get rid of doubling 
contribution of the fermionic propagator. 
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